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We present an ab initio analysis of electron conduction 
through a Ceo molecular device. Charge transfer from the 
device electrodes to the molecular region is found to play 
a crucial role in aligning the lowest unoccupied molecular 
orbital (LUMO) of the Ceo to the Fermi level of the elec- 
trodes. This alignment induces a substantial device conduc- 
tance of ~ 2.2 x (2e 2 jh). A gate potential can inhibit charge 
transfer and introduce a conductance gap near Ef, chang- 
ing the current-voltage characteristics from metallic to semi- 
conducting, thereby producing a field effect molecular current 
switch. 
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Understanding electron conduction in atomic and 
molecular scale nano-devices is an extremely active-re.- 
search topic at present. In several recent experimentsETH, 
electric current flowing through truly atomic and molec- 
ular scale conductors, as small as a single atom, have 
been measured. The current voltage (I-V) characteris- 
tics of these molecular devices have profound potential 
for device application, including high nonlinearity, neg- 
ative differential resistance, and electro-mechanic cur- 
rent switching. However, a thorough understanding of 
electron transport mechanisms at this scale, which will 
play a crucial role in designing operation principles for 
future nanoelectronics, has not yet been achieved. To 
provide insight into electron conduction mechanisms at 
the atomic and molecular scale, we have investigated the 
quantum transport properties of a Cqq based molecular 
device using a first principles theoretical method. 

Experimentally, the Cqq molecular electro-mechanical 
amplifier^ demonstrated that electric current flowing 
through a Cqq molecule could be amplified by as large as 
100 times when the molecule was mechanically deformed. 
This phenomenon was attributed to an increase of the 
tunneling density of states at ihe Fermi level induced 
by the mechanical deformatioroLl. While this conduc- 
tion mechanism is very interesting, a natural question is 
whether or not there are other, perhaps simpler, physical 
mechanisms which can be exploited for molecular scale 
electron conduction ? Our analysis suggests a unique and 
new conduction mechanism based on charge transfer dop- 
ing (CTD): because of its high electro-negativity, a well 
contacted Ceo gains charge from CTD thereby aligning 
its LUMO state to the Fermi level of the electrodes re- 
sulting in a substantial device conductance without the 
need of mechanical deformation. CTD was seen to have 



reduced the conductance of a short carbon chainla, but 
for the more complex Cqq molecule, it does the opposite 
by drastically increasing conduction. Current switching 
mechanisms are another important yet unsettled ques- 
tion of nanoelectronics. Our analysis further suggests 
a new principle for molecular switching: using proper 
electro-negative molecules and metallic electrodes, cur- 
rent can be switched on and off by controlling the amount 
of charge that is transferred to the molecular region. 

The 6*60 device we study is illustrated schematically 
in Fig. (la) where a Cqo molecule is electrically well con- 
tacted by two atomic scale Al metallic electrodes which 
extend to reservoirs far away where bias voltages Vi >r are 
applied to the left (I) and right (r) electrodes, respec- 
tively. An additional gate voltage V g may also be applied 
to a metallic gate capacitively coupled to the molecule. 
Due to the large number of atoms and complications 
such as the gate, the bias, the atomic electrodes (as op- 
posed to jellium electrodes), and the presence of localized 
states, existing ab initio methods for analyzing quantum 
transporttTO cannot be applied. We have therefore de- 
veloped a new approach which combines non-equilibrium 
Green's function theorytjO with pseudo-potential real 
space ab initio simulation techniqueaU. It is worth not- 
ing that conventional density functional theory (DFT) 
methods solve problems for cither finite systems such as 
an isolated molecule, or periodic systems consisting of 
super-cells. In contrast, central to our problem of pre- 
dicting quantum transport properties is the ability to 
deal with systems having open boundaries provided by 
long electrodes which maintain different chemical poten- 
tial due to external bias. In other words, a typical device 
geometry, e.g. Fig. (la), is neither isolated nor periodic. 
This situation is most convenie ntly handled by nonequi- 
librium Green's function theoryOO. ._. 

Briefly, our technique is outlined as follows _ . We di- 
vide the long device system into three regionst2l: the left 
and right electrodes, and the scattering region (Fig. (la)). 
The scattering region actually includes a portion of 
the semi-infinite electrodesEHl. After self-consistency is 
reached, this arrangement ensures good electric contact 
between the electrodes and the molecule, establishes a 
common Fermi level for the system, and ensures charge 
neutrality at .equilibrium. When the scattering region is 
large enoughEJES, the Kohn-Sham (KS) potential out- 
side the scattering region is well screened and therefore 
well approximated by a perfect "bulk" electrode epyiron- 
ment which we obtain by a separate calculation!!]. By 
setting the KS potential outside the scattering region to 
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the bulk value and matching the electrostatic potential 
at the boundary (including the gate), the infinite open 
boundary problem is reduced to the calculation of charge 
density inside the finite scattering region with the elec- 
trodes' contribution accounted for by appropriate self- 
energiesy-OEI 

To perform our self-consistent DFT analysis, we calcu- 
late the charge density from the density mat rix w hich is 



related to nonequilibrium Green's functio 
tor, 
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and fi ir {E\ in >r ) are the distribution functions deep in 
the left and right electrodes with chemical potentials 
Hi, r , respectively. The reservoirs are assumed to be at 
equilibrium, therefore fi tT ks 0(E — /i;. r ) at low tem- 
peratures where O is the step function. The retarded 
(and advanced) Green's functions G R (G^V-ofj-the sys- 
tem are then calculated in standard fashior£lTt£l by di- 
rect matrix inversion where self-energies due to coupling 
to electrodes (Sfi, S^ r ) are calculated by extending the 
method of Ref. [18| Because there are many bands iftpur 
atomic electrodes (as compared to Jellium electrodesoEj) , 
the energy integration of Eq. (fy) must be done with 
care so that global charge neutrality is maintained at 
equilibrium. With the charge density, we can evalu- 
ate the effective device potential V^g[p(r)] which consists 
of Hartree, exchange-corxelation, atomic core, and any 
other external potentialsEil. Using the boundary condi- 
tion on the effective potential discussed above, V c g(r) 
is defined everywhere in space. It is therefore straight- 
forward to construct the Hamiltonian-jnatrix using an 
s, p Fireball atomic orbital basis setEJ. We use stan- 
dard norm-ceaserving pseduopotcntials to describe the 
atomic coresE3. The DFT iteration is repeated until 
self-consistency, and current is c alculate d by integrat- 
ing conductance G(E) over energyliSEjO, with G(E) = 

AG Tr[/m(Sf )G fl /m(Sf r )G A ] where G a = 

Fig. (lb) shows the equilibrium charge_.distribution 
along the middle cross-section of the deviceO at zero bias 
and gate voltages. The charge density in the electrodes 
is affected by the 6*60 but this effect is well screened 
away from the molecule. Indeed, the charge density con- 
tours match perfectly at the connections between the 
electrodes and the scattering regionEJ. The electronic 
structure of an isolated Cqq- molecule is well knownc3 and 
our ab initio calculation gives a HOMO-LUMO gappof 
1.77eV, in excellent agreement with previous literature^. 
The undoped Cgo molecular solid is a semiconductor£3. 
Since Cqq is electro-negative, one can dope the solid and 



fill the LUMO state with up to 6 electronso. Impor- 
tantly, for our device the open metallic electrodes pro- 
vide natural doping through CTD. As a result, we found 
that, in equilibrium, three extra electrons flow into the 
molecular region. This is a substantial charge transfer in 
order for the molecule to equilibrate with the electrodes 
so that a common Fermi level is established. Experi- 
mentally, it is found that the Cgo solid conducts best 
when it is doped with three electrons per CqoE3. This 
is because doping three electrons amounts to half filling 
the LUMO. Our ab initio analysis predicts, for the Cqq 
device, that the Cqq LUMO state is half-filled and we 
find an equilibrium conductance of G(Ep) = 2.2G at 
zero temperature. This value is very significant: without 
charge transfer we would expect a much smaller conduc- 
tance due to the filled HOMO state and the substantial 
HOMO-LUMO gap of the isolated C 60 . 

In Fig. (2), we plot the equilibrium transmission eigen- 
values of the device as a function of electron energy at 
two different gate voltages. At V g = 0, there are three 
transmission eigenvectors which contribute substantially 
at Ep (the vertical line) to the equilibrium conductance 
G = 2.2G . However, a gate voltage shifts the states 
near Ep and changes the transmission channels signifi- 
cantly. This is shown in Fig. (2b) where we see a conduc- 
tance gap at Ep. Therefore, a gate potential can change 
electron conduction through this device significantly, pro- 
ducing field induced molecular switching. We can under- 
stand the transport properties by projecting the scatter- 
ing states of the device on to the molecular orbitals of an 
isolated Cqq- We found that scattering states near Ep 
are over 90% HOMO or LUMO in character, indicating 
that the Cqo retains its isolated electronic structure and 
that other Cqo orbitals are less important for conduction. 
In addition, we can classify each scattering state as being 
majority HOMO or majority LUMO. We find that ma- 
jority HOMO states are almost pure HOMO (> 90%); 
while majority LUMO states are a HOMO-LUMO mix- 
ture (~30:70). Specifically, the transmission eigenvec- 
tors of Fig. (2a) are all LUMO-like, indicating that there 
are three electrons half filling the LUMO state and con- 
tributing three conduction channels at equilibrium. In 
other words, the LUMO state is naturally aligned with 
Ep due to charge transfer. On the other hand, a negative 
gate voltage can inhibit charge transfer so that Ep lies 
between the HOMO and LUMO states of the molecule. 
Indeed, we found that the scattering state corresponding 
to the left peak in Fig. (2b) is HOMO-like, and that of the 
right peak is LUMO-like. In the absence of charge trans- 
fer, one would expect transmission properties similar to 
those in Fig. (2b). 

The above physical picture is further demonstrated in 
Fig. (3) which shows the equilibrium conductance G, the 
three individual transmission eigenvalues Ti (i = 1,2,3), 
and the number of transferred charge Q (inset) as func- 
tions of gate voltage V g at Ep. G has a step-like behav- 
ior: each "step" is due to the depiction of a LUMO-like 
state. Only after charge transfer is completely inhibited 
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do we obtain the conductance gap illustrated in Fig. (2b) 
where G becomes very small. Our result therefore indi- 
cates that charge transfer plays a crucial role and changes 
the physical picture of electron conduction in this device 
qualitatively. The behavior of G also has a clear corre- 
spondence with Ti (the dotted lines with circles): as V g 
is scanned toward more negative values, TVs are removed 
one by one. 

Charge transfer also has very important implications 
for the I-V curves, shown in Fig. (4). The current is plot- 
ted as a function of right bias V r fixing Vi = 0. When 
V g = the I-V curve shows clear metallic behavior. 
When V g ^ 0, the Ceo device can change from a metal 
(Fig. (2a)) to a semiconductor (Fig. (2b)) with a con- 
ductance gap near zero bias. Therefore the gate potential 
can "switch" off the current, reflecting the transition of 
transmission eigenvalues from Fig. (2a) to (2b). This cor- 
responds to an amplification factor of ~ 20. The value of 
gate potential is non-universal as it depends on the shape 
of the gate and dielectric medium surrounding the device. 
What is essential is to produce enough electric field lines 
inside the molecular region so that the electrode-doped 
LUMO electrons are depleted from the molecular junc- 
tion. For our system, a shift of ~ O.leV is enough to 
generate switching. This is however a large shift and it 
suggests that a sharp-shaped gate is perhaps necessary 
in an experimental setup, otherwise most field lines will 
be screened by the metallic electrodes. The predicted 
current is in the range of ~ 10fj,A at small bias volt- 
ages (~ 50mV). This is somewhat larger thaflr-the ~ 4/xA 
measured for the electro-mechanic amplifierEJEJ. 

Our results also have important implications and pro- 
vide a benchmark for semi-empirical and other theo- 
ries, such as those based on parameterized tight bind- 
ing models. The chemical potential difference between 
an isolated Al electrode and an isolated Cqo molecule 
is A = Ef ectrode - E% 60 ~ O.Wa.u.. This provides a 
lot of uncertainty as how one should align the molec- 
ular orbitals to the Fermi level of the electrodes in a 
non-self-consistent calculation. SpecificallyO, calculat- 
ing Q(^Ep ectrode ) without alignment gives 1.17G , much 
less than the correct result. Shifting levels of Cqq by 
adjusting a chemical potential so that global charge neu- 
trality is obtained, we get 1.7G . Shifting levels so that 
three extra charges are inside Cqq region, the result be- 
comes 2.0G o . Because we cannot know charge transfer 
in a non-self consistent analysis, the result suggests that 
requiring global charge neutrality in a semi-empirical cal- 
culation is the next best thing to do at equilibrium. In 
this regard, we note that there has been some practice 
for requiring local charge neutrality on each atomcScIl 
in semi-empirical calculations. Another interesting re- 
sult we found is that there are many localized states 
in the molecular junction. At equilibrium, we find 96 
bound states by integrating the density of states (Eq.(jl])) 
from — oo to the propagating threshold of the electrodes. 
These localized states play an important role in estab- 
lishing the effective potential V c s and must be included 



in the analysis. 

In summary, our results suggest that charge transfer 
doping induces a substantial conductance in a well con- 
tacted molecular Cqo device. Essentially, CTD aligns 
the LUMO states of the Cqo with the Fermi level of 
the electrodes, opening up 3 conductance channels and 
producing metallic I-V characteristics at equilibrium. A 
field effect provided by a gate potential can switch off 
the current by inhibiting charge transfer. This is very 
interesting from a device operation point of view and 
provides a new mechanism for molecular scale current 
switches, complementing the electro-mechanical opera- 
tion principle explored in previous worka. The charge 
transfer induced conduction discussed here should be a 
more general mechanism: using proper electrodes which 
provide natural doping to electro-negative molecular sys- 
tems with a filled HOMO state and substantial HOMO- 
LUMO gap, an increase in conduction is expected. Since 
one expects charge transfer to be a generic feature of 
molecular devices, this effect should be exploited further 
for applications of molecular scale current switches. 
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FIG. 3. Solid squares: equilibrium conductance G(Ef) 
as a function of gate voltage V g . Vertical dotted lines indi- 
cate, approximately, integer number of transferred charges Q. 
Open circles: the three transmission eigenvalues as a function 
of V g at Ef- Inset: the transferred charge as a function of 
gate voltage. 

FIG. 4. I-V curves of the system showing the switching 
between metallic (squares) and "semiconducting" (circles) be- 
havior. 



FIG. 1. (a). Schematic plot of the Cqo molecular device, 
(b). Contour plot of the equilibrium charge density. Notice 
the perfect match across the boundaries between the scatter- 
ing region and the electrodes. 



FIG. 2. Transmission eigenvalues as a function of electron 
energy, (a) For V g = where there are three transmission 
eigenvalues, (b) For V g = la.u. Vertical line shows the Fermi 
level of the system. 
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